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ABSTRACT 



The BATSE mission aboard CGRO can be used to 
observe hard X-ray sources by using the Earth oc- 
cultation method. This method relics on measuring 
a step in the count rate profile in each BATSE detec- 
tor as a source rises above or sets below the Earth's 
limb. A major problem in determining the step sizes 
(and hence the flux) is in extracting the steps from 
the varying background. A technique for flat-fielding 
the response of 7-ray detectors has been developed at 
Southampton. The technique uses a dynamic Monte- 
Carlo model to simulate the dominant components 
of the gamma-ray background encountered by the 
experiment at any point in its orbit. A maximum 
likelihood imaging method is also being developed 
that will be used to make a sky survey with all ~9 
years of the BATSE CONTINUOUS data set in the 
20 - 500 keV range. An all sky map of 25 - 35 keV 
emission has been made using 60 days of data and 
has a 3cr flux sensitivity of 14 mCrab. 
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1. INTRODUCTION 



The Burst and Transient Source Experiment 
(BATSE) consisted of eight 2025 cm 2 Nal(Tl) Large 
Area scintillation Detectors (LADs), sensitive from 
20 keV - 2 MeV. BATSE monitored the whole sky 
continuously from April 1991 - June 2000 as an ex- 
periment aboard the Compton Gamma-Ray Obser- 
vatory (CGRO). The continuous whole sky exposure 
makes BATSE an ideal facility for the production of 
the first all sk y hard X-ray survey since HEAO A-4 
of 1977-1978 ( |Lcvinc ct al. (1984J )). A survey per- 
formed with BATSE will not only be the first in the 
few hundreds of keV band but will have the ability 
to detect weak persistent sources and continuously 
monitor variable and transient sources. The BATSE 
CONTINUOUS data set contains usable count rates 
in 15 energy channels recorded in 2.048 s time bins. 



BATSE EARTH OCCULTATION 
OBSERVATIONS 



BATSE was developed primarily to detect and locate 
7-ray bursts. The individual BATSE detectors have 
no positional sensitivity hence to locate the position 
of a source of 7-rays the relative count rates from 
several detectors are combined, alo ng with the point 



spread function of each detector ( Pendleton et al 
(1999| )). Prior to launch however it was realised that 



that the instrument could also be used as a sensitive 
all sky monitor of point like 7-ray sources by using 
the Earth to modulate the flux at the CGRO orbital 



period ( Paciesas et al. (1985 )). Sharp steps are seen 
in the count rate profile of the detectors as point 
sources rise and set over the Earth's limb, the si ze 
of which give the source flux ( |Harmon et al. (2000| )). 
About 70% of the whole sky is occulted twice per 
90 minute orbit. The precession of the orbit is such 
that the whole sky undergoes occultation every 53 
days. The key to measuring the flux values correctly 
is in accura tely determining the varying background 
count rate. Ling et al. (200C ) have recently published 
a catalogue of 34 moderately strong 7-ray sources 
observed with BATSE during the first three phases 
of the CGRO mission (May 1991 through October 
1994). They use an empirical fit to find the back- 
ground, which is hamper ed by systematic erro rs and 
large computation time ( Harmon ct al. (2000] )), and 
make no attempt to image the sky. 



SOUTHAMPTON BATSE BACKGROUND 
MODEL 



The BATSE model is based on The INTEGRAL 
Mass Model (TIMM), which is a GEANT based 
Monte-Carlo simulation code devloped at Southamp- 
ton to investigate the performance of the instru- 
ments aboard th e INTEGRAL mission pre-launch 
( Lei et al. (1999 )). It considers three components to 
the background; cosmic rays, cosmic diffuse gamma 
rays and 7-rays from interactions of cosmic rays in 
the Earth's atmosphere. These same techniques are 
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Figure 1. Comparison of the count rate profiles recorded by the BATSE LAD 5 (dots) and the Southampton 
modelled data. The total is the sum of the cosmic diffuse emission (CD), the atmospheric albedo (ATM) and 
prompt cosmic ray interactions (CR). 



easily extended to other experiments and has been 
proposed for use with the BAT burst monitor aboard 
SWIFT. 



3.1. Dynamic Generation of Background 
Components 



With BATSE we are interested in the dynamical 
behaviour of the background within the CGRO or- 
bit. The background is a complicated function of the 
pointing of the detectors, the position of CGRO in 
the Earth's magnetic field and its irradiation history. 
A novel two step method has been developed to ac- 
count for the background at any point in the CGRO 
lifetime. Firstly CGRO is assumed to be in deep 
space without any interference from Earth. Here it 
receives isotropic fluxes of cosmic rays, cosmic diffuse 
7-rays and atmospheric 7-rays and a database, con- 
taining information about all the events interacting 
in the BATSE modules is created. The second stage 
filters the database with the physical criterea cor- 
responding to CGRO's position and pointing. This 
information is used to create the background for any 
orbital position and pointing. The component due to 
the atmospheric 7-rays is found by scaling and sam- 
pling of the detector pointings and finally the radio- 
activation of the detector materials is simulated. The 
first stage is computationally very intensive, but only 
needs to be done once for the whole lifetime of the ex- 
periment. To compute the second stage of the back- 
ground simulation requires only a few hours for all 
of the eight BATSE modules for one day. Figure |l| 
shows the excellent agreement between the modelled 
data and the count rate profile in the BATSE LADs 
for three energy channels. 



Table 1. Approximate sensitivity and angular res- 
olution of the imaging technique at ~30 keV as a 
function of filter width, f . 



f 


3fi flux 


Angular resolution 


(bins) 


sensitivity (mCrab) 


(HWHM°) 


5 


270 


0.5 


20 


199 


0.9 


40 


88 


2.4 


GO 


77 


2.9 


120 


74 


5.3 



4. PRODUCTION OF SKY IMAGES 



Sky images are produced using a maximum likeli- 
hood meth od originally develo ped by J. Knodlseder 
at CESR (iKnodlseder (1999| )). This method has 
been adapted at Southampton for use with the 
flat-fielded CONTINUOUS data and developed into 
the LIMBO code (Likelihood Imaging Method for 
BATSE Occultation). 



4.1. Flat-Fielding of Data 

The data is flat-fielded by subtracting the modelled 
background and the result passed to a simple differ- 
ential filter of width / bins, given by 



d' = 



(i) 



The filtered data, d[, is transformed from a count 
rate containing occultation steps to a set of occul- 
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Figure 2. Map of the 25 - 35 keV flux from the whole sky for a total of 61 days from TJD 10840 to TJD 10900 
made using LIMBO. The colour scale of the map has been limited to 250 mCrab (one quarter of the peak of 1 
Crab) to show details at lower significance. Crosses mark the actual positions of the objects given in table^ 



tation peaks. A square step at bin i will be trans- 
fromed to a triangular peak of width 2f. The choice 
of / has implications for the sensitivity and angular 
resolution of the imaging technique; larger / will give 
a more significant occultation peak and improve sen- 
sitivity, smaller / returns less significant but sharper 
peaks that will lead to better angular resolution. The 
gain in sensitivity from increasing / begins to level 
out as the steps from other sources enter the filter. 
Table @ summarises the sensitivity and angular res- 
olution obtained from images made from one day's 
observation of Cygnus X-l using energy channel 2 
(^25 - 35 keV). The creation of images is discussed 
further in section [4.2| . The raw B ATSE data contains 
occasional gaps where, for example, the instrument 
is turned off during passage through the South At- 
lantic Anomaly (SAA) or where contact is lost be- 
tween CGRO and telemetry relay satellites. Left un- 
treated the gaps cause large artificial peaks in the 
data near each gap so a linear interpolation is used 
to fill the gaps in the background subtracted data. 
The gap data is not used further in the analysis. 



4.2. Maximum Likelihood Imaging 



Zhang et al. (1995 )) in that it produces the results of 
statistical tests for emission at each point in the map 
rather than a deconvolution of the data. The occul- 
tation peak dataset, o, is fitted to a response vec- 
tor, e, which describes the expected position, shape 
and amplitude of the peaks as a function of source 
strength and position on the sky. The fit determines 
a scaling factor of the response vector, which is pro- 
portional to the flux received from each point in the 
sky for all eight LADs. The scaling factor, a, is given 

by 



a 



(2) 



where Cj is the statistical error on Oj. A maximum 
likelihood ratio test is used to determine the signif- 
icance of the flux detection and ret urns a skyma p 
for a given grid of source positions (Cash (1976)). 
The maximum likelihood ratio, A, is distributed as 
xf, with v degrees of freedom. It is calculated from 
the difference C — C src , the \ 2 statistics for the null 
(no source is fitted to o) and source (a source with 
free strength is fitted to o) hypotheses respectively, 
where 



The imaging method used here differs from the 
Radon transform method used by other groups (e.g. 



C 



E 



(3) 
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(oi - aei) 2 



(4) 



When searching for a known source v = 1 and hence 
the detection significance, in gaussian er, is given by 
A /A. For searches of new sources of 7-rays v = 3 



(Knodlseder (1999)). The image is built up from the 
superposition of the arcs caused by the occultation 
of point sources by the Earth's limb. The angular 
resolution perpendicular to the limb depends on [3, 
the angle at the geocentre between the orbital plane 
and the source. Sources seen with (3 ~45° yield good 
information since the rising and setting limbs are al- 
most orthogonal. 



The maximum likelihood method has an advantage 
over the Radon transform method in that there is 
no limit to the size of map that can be generated. 
The Radon transform is not well behaved for regio ns 
much larger than ~20°x20° ( [Harmon ct al. (2000[ )). 



Likelihood maps can be made for short periods of 
time and then added together to increase source sig- 
nificance whilst retaining the ability to investigate 
variability on short timescales. The total likelihood 
ratio from a total of individual daily maps, indexed 
by d, is A T = J2d(^d - A d ) where 

A d = (o T - a d ){{a T - a d )D d - 2F d ) (5) 

and F d and D d are the daily values of the numera- 
tor and denominator of equation respectively. The 
quantity olt is the total flux for all days being con- 
sidered, i.e. a T = J2d a d = J2 d F d/J2d D d- 



Table 2. A selection of point sources of 25 - 35 keVj- 
rays found from the sky map of figure with flux > 50 
mCrab. Sources were included only if the distance 
between the peak flux and the actual source position 
(final column) was < 3°. 



Source 


Peak 


Flux 


Peak 


Distance 


Name 


a 


(mCrab) 


(tf>) 


(degs) 


Crab 


212 


1000 


( -177, -6) 


1.58 


CygX-1 


154 


727 


( 71, 4) 


0.99 


ScoX-1 


71 


338 


( -1, 24) 


0.24 


4U1700 


G7 


319 


( -13, 2) 


0.77 


1E1740 


G2 


294 


( -1, o) 


1.05 


OA01657 


60 


282 


( -15, 2) 


1.44 


GX354+0 


58 


273 


( "5, 0) 


0.72 


VelaX-1 


49 


231 


( -99, 4) 


2.06 


4U1630-47 


40 


188 


( -25, 0) 


1.93 


CygX-3 


38 


181 


( 77, 0) 


2.92 


GRS1915 


28 


132 


( 43, 0) 


1.94 


AqlX-1 


20 


94 


( 35, -2) 


2.26 


GX301-2 


20 


93 


( -61, 0) 


1.10 


Mrk501 


18 


84 


( 61, 40) 


2.82 


PKS2005 


17 


80 


( 11, -30) 


0.81 


Geminga 


16 


77 


( -163, 2) 


2.94 


4U1608 


16 


77 


( -29, -2) 


1.15 


CenA 


14 


G7 


( -49, 18) 


2.06 


3C273 


16 


74 


( -69, 62) 


2.58 


HR1099 


14 


65 


( -175, -40) 


1.57 



5. RESULTS 



Figure |2| shows the total likelihood map for the whole 
sky created from channel 2 CONTINUOUS data for 
60 days between Truncated Julian Days (TJD) 10840 
- 10900. The energy calibration of BATSE channel 2 
is such that the range of photon energies it detects is 
approximately constant at ~9 keV with a lower limit 
of ^25±5 keV. The Crab, Cygnus X-l and several 
Galactic Centre sources have been identified along 
with other faint sources as summarized in table 0. 
Some artifacts remain around the brighter sources, 
which we aim to reduce by removing the brighter 
sources iteratively. The 3cr flux sensitivity of the 
map is 14 mCrab. The map was created using a 
2°x2° sky grid and / = 20. Table | suggests that 
it would be worthwhile increasing the grid resolution 
to 0.5° xO. 5°, which would increase the computation 
time and disk space requirements by a factor of 8. 
The current maps are ~600 kB in size per day and 
energy channel and LIMBO takes ~45 mins to run on 
a Pentiumlll Linux workstation, although the code 
is not fully optimised as yet. On the basis of this 
preliminary analysis LIMBO shows a lot of promise 
as a method of surveying the sky. We will continue to 
develop LIMBO with the aim of creating a database 
of daily maps in each energy channel for the entire 
BATSE observation. Sources found in the surveys 
will be studied in more detail for location, variability 



and spectroscopy (see also presentations by Dean and 
Westmore at this meeting). 
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